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We consider the classical 0(jV)-symmetric models confined in a d-dimensional slab-like geometry 
and subject to periodic boundary conditions. Applying the one-particle-irreducible variant of func- 
tional renormalization group (RG) we compute the thermal Casimir forces acting between the slab 
boundaries at criticality. The applied truncation of the exact functional RG flow equation (leading 
order in the derivative expansion) neglects renormalization of the momentum dependencies in the 
propagator, but retains interaction vertices of arbitrary order. We evaluate the critical Casimir 
amplitudes A/(d, N) for continuously varying dimensionality between two and three and AT = 1, 2. 
Our findings are in very good agreement with exact results for d = 2 and N = 1. For d = 3 our 
results are closer to Monte Carlo predictions than earlier field-theoretic RG calculations. 
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I. INTRODUCTION 

The concept of Casimir forces^- i.e. long-ranged 
effective interactions between macroscopic bodies im- 
mersed in strongly fluctuating media, is nowadays rec- 
ognized in plethora of systems ranging from biology to 
cosmology. On one hand their prominent importance 
stems from providing possibilities of testing theories of 
fundamental interactions^— on the other they are in- 
voked to explain important mechanisms underlying phys- 
ical phenomena in organic matter—^ and as a con- 
trol tool for nanodevices^ The more conventional and 
well-studied (also experimentally) examples of Casimir 
forces in condensed matter include fluid mixtures^ 3 - - — 
and helium J£r— 

The essence of the phenomenology underlying the 
Casimir forces is the presence of soft fluctuations occur- 
ring within a medium, whose spectrum of excitations be- 
comes constrained by the boundary conditions imposed 
by the confining walls or by macroscopic objects im- 
mersed therein. As a result, the free energy acquires a 
contribution depending on the separation L between the 
walls (or the objects) and the system finds it favorable 
to either increase or decrease this distance. A promi- 
nent characteristic of Casimir forces is the universality 
of its asymptotic behavior (including the amplitude) for 
L ^s> Aq" 1 and £ 3> L, where A^" 1 is a system-specific 
microscopic lengthscale and £ denotes the (bulk) corre- 
lation length for the medium fluctuations. On the other 
hand, the amplitude of the force (even its sign) does de- 
pend on the boundary conditions. The orthodox QED 
Casimir forced is a pure effect of quantum vacuum fluc- 
tuations. Its theoretical description involves a free field 
theory, where interactions with matter are taken into ac- 
count via boundary conditions and the expectation value 
of the fluctuating electromagnetic field is zero. The situa- 
tion is quite distinct in these respects in many condensed 
matter systems. Here a field-theoretic description of the 
problem often requires a treatment of interaction terms 
in the effective action. For certain cases, the boundary 



terms may also lead to nonzero expectation values of the 
(nonuniform) order parameter, which gives rise to effec- 
tive interactions between the confining walls (or other 
macroscopic objects) also in the absence of fluctuations. 

In this work we analyze the thermal Casimir effect em- 
ploying a class of standard models to describe different 
aspects of critical behavior, namely the d-dimensional 
0(./V)-symmetric ^ 4 -type models. The system is con- 
fined to be finite in one direction (so that its extension 
is L) and we consider the periodic boundary conditions. 
Despite not being of the highest relevance for quanti- 
tative comparison with experimental data, this choice 
of boundary conditions does not yield Casimir forces at 
mean-field level, so that the effective interaction is a pure 
fluctuation effect (preserving the analogy to the QED 
setup). In addition, these boundary conditions yield a 
zero-mode in the free propagator, which leads to substan- 
tial complications in a conventional RG treatment, where 
the e-expansion2£ becomes ill-defined beyond two-loop 
orde n 26 i 27 We study the system within the one-particle- 
irreducible (Wetterich) framework of functional RG. 28 i 29 
The exact flow equation governing the flow of the gener- 
ating functional for irreducible vertices is treated within 
an approximation retaining the flow of interaction cou- 
plings of arbitrary order, but neglecting the renormal- 
ization of the momentum dependencies in the propaga- 
tors. This amounts to the leading order of an approxima- 
tion strategy known as derivative expansion, 2 - which was 
successfully applied inter allia for accurate calculations 
of the O(N) bulk critical behavior— ~— (including the 
Kosterlitz-Thouless case22). This approach is not con- 
trolled by any small parameter, but may be considered 
as a leading order of a systematic procedure to compute 
critical properties. It can be applied for different d and 
N. In particular, we scan the dependence of the critical 
Casimir force amplitude Af(d,N) as a function of di- 
mensionality d, varying d between two and three at fixed 
N G {1,2}. Our findings compare very well to the exact 
valued for d = 2 and N — 1 and reasonably well with 
numerical Monte Carlo (MC) simulation s 34 ' 35 in d = 3 
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(both for JV = 1 and N — 2). 

The paper is organized as follows: In Sec. 2 we present 
the model and the RG setup. In Sec. 3 we consider the 
bulk critical behavior upon varying d between two and 
three at fixed number of order-parameter components 
JV G {1,2}. In particular, within the applied approxima- 
tion, we compute the correlation length in the vicinity 
of the critical point and extract the values of the critical 
exponent v{d, JV), describing the divergence of the corre- 
lation length £ at criticality. In the subsequent Sec. 4 we 
present the computation of the critical Casimir forces in 
the cases JV = 1 and JV — 2, continuously varying d. We 
summarize the work in Sec. 5. 



II. MODEL AND RG SETUP 

The Landau-type model considered here is a paradigm 
of the theory of phase transitions and critical phenom- 
ena. In particular, in Refs. I25U27I it was employed to 
study the thermal Casimir forces. The fluctuating N- 
component real order-parameter field x(x) i s confined in 
a d-dimensional box of volume V — D d ~ x x L, where D 
is much larger than all the other lengthscales present in 
the system. We impose periodic boundary conditions in 
all the directions and the effective Hamiltonian takes the 
standard form 

S{x\ = J d d x[U Q ( P ) + ^o(Vx) 2 ] , (1) 

where we introduced p — \{x) 2 ■, and the effective poten- 
tial 



U (p) = ^ + S P 



(2) 



is quadratic in p. The model involves the interaction 
coupling v, the mass-like parameter 8, the constant Zq, 
and is additionally supplemented with a lower cutoff in 
real space at Aq , which we treat on equal footing with 
the other parameters. 

The present approac h 29 i 36 relies on an exact RG equa- 
tion governing the flow of the scale-dependent effective 
action I\ [</>], where the field </> is the expectation value 
of x(x)£2- Upon decreasing the momentum cutoff scale 
A from Ao to zero, the functional Ta^] interpolates be- 
tween the initial action S[x — <fr] a ^ A — Ao and the 
full effective action, i.e. the free energy To[4>] at A = 0. 
The quantity T\[(j)] may be understood^ as the free en- 
ergy obtained after including the fluctuation modes with 
momentum q G [A, Ao]. At A = Ao no fluctuations are 
included and T &[(/>] is just the bare effective action. On 
the other hand, in the limit A — > all the fluctuation 
modes are included. The flow of Ta^] upon reducing 
the cutoff scale A is governed by the flow equation^ 



= - 2 ^ m 



(3) 



where R A (q 2 ) is a momentum cutoff function added to 
the inverse propagator to regularize the infrared behav- 
ior, and r^[$] = 5 2 T A [$\/ {5^8^), with a G {1, ...JV} 
(no off-diagonal components of the inverse propagator oc- 
cur in the present problem) . The trace Tr sums over field 
components and momenta: Tr = ^~J Q j- ^2 qi J (2^)^-1 m 
the limit D — > 00 and with qi = n G Z . Inte- 

gration of Eq. (lU) with the initial condition at Ao spec- 
ified by Eqs. (TTJ) and © yields the exact solution to 
the model. Herein we apply an approximation strategy 
known as derivative expansion,— where the symmetry- 
allowed terms in Ta [4>] are classified according to increas- 
ing number of derivatives and we keep only the leading 
(zeroth) order term in this expansion. This truncation 
level (the so-called local potential approximation (LPA)) 
amounts to imposing the following ansatz: 



d d x[U A (p) + \z (v$f 



(4) 



where we put p = ^(4>) 2 . The present approximation 
therefore neglects renormalization of all the momentum 
dependencies in the propagator, but does not truncate 
interaction terms in the local effective potential U\(p) at 
any order. 

With this ansatz the effective action flow is represented 
by the flow of the scale-dependent effective potential 



d A U A (p) = ^Tr \ d A i? A (<T 



f 



JV- 1 



M R (p,q 2 ) M T {p,q 2 ) 



where 

M R (p,q 2 ) 
M T {p,q 2 ) 



(5) 



Z q 2 + R A (q 2 ) + U' A (p) + 2pU'l{p) , 
Z q 2 + R A (q 2 ) + U' A (p) , (6) 



and Tr = £ J S^yA ■ Since we neglect the renor- 
malization of Zq, Eq. ([5]) is closed and our subsequent 
analysis is based on its numerical integration. Any prac- 
tical calculation requires specifying the cutoff function 
R A (q 2 ), which must respect a number of conditions^ 
We make the following choice, — known as the Litim cut- 
off 



R A {q 2 ) = Z (A 2 - q 2 )6(A 2 - q 2 ) 



(7) 



for which the trace in Eq. ([5]) can be straightforwardly 
performed in the limit L — > 00. 
The free energy is identified as 

F{<f>,L) = fc s rr A ^ o [0 = <Punif) = k B TVU A ^ (p) , (8) 

where the functional is evaluated at a uniform field con- 
figuration. Its dependence on L arises due to the trace in- 
volving the L-dcpendent momenta allowed by the bound- 
ary conditions. The Casimir forces, in turn, result from 
the L-dependent free energy. Since our study concerns 
the long-ranged interactions close to a second-order phase 
transition, we first need to resolve the critical behavior, 
which is done in the next section. 
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III. CRITICAL BEHAVIOR 

In order to efficiently tune the system to its bulk crit- 
icality, we consider the limit L —> oo and perform the 
standard variable transformation 



p = ZoA 2 - 1 



P 



u%(p)=A- d U A (p) , 



(9) 



where the superscript oo indicates that we deal with an 
infinite system. This brings Eq. ([5]) to a scale-invariant 
form 



1 N -1 



A, 



mn.{p) mr(p) 



(10) 



where 



m R (p) = l + uT'(p) + 2puT"(p), 



m T (p) = 1 



'00- 



(11) 



The primes denote differentiation with respect to the 
argument. We introduced s = — log(A/A ) and Ad 

where S d is the surface area of a d-dimensional 



(2TT) d d ' 

unit sphere. 

We numerically integrate the discretized version of 
Eq. (fTUj) taking Eq. (T2]) as the initial condition. For all of 
the analysis we put Zq = 1 and Ao = 1. At fixed d and 
the (initial) interaction coupling v in Eq. (T2|), we scan the 
solution to Eq. (fTt)|) for A — > upon varying the bare mass 
parameter 6 in Eq. (f2|). We identify two regimes in the 
space spanned by (5, u), corresponding to the symmet- 
ric and the symmetry-broken phases. The line of critical 
points (in the (6, v) space) and the corresponding fixed 
points are found by the dichotomy procedure illustrated 
in Fig. 1. Having identified the critical point, we vary <5 
so that the system remains in the symmetric phase close 
to S c and scan the behavior of the (non-rescaled) effective 
potential at its minimum in the limit A — > 0. By com- 
puting the curvature of the effective potential as function 
of <fi at <fi — for different values of 6 > 8 C in the vicin- 
ity of the critical point, we extract the correlation length 
£(6). The obtained behavior very well fits the expected 
power- law 



(12) 



We analyze the dependence of the exponent v(d, N) fo- 
cusing on N e {1,2} and varying d between two and 
three. The obtained dependencies are shown in Figs 2 
and 3 for N = 1 and N = 2, respectively. In both 
cases the values of v compare well to the best available 
estimates^ in d = 3. We obtain 2/(3, 1) « 0.64 and 
2/(3, 2) ss 0.68 (the best estimates according to Ref. [H 
are 2/(3, 1) ~ 0.63 and 2/(3,2) » 0.67 up to two digits). 
On the other hand, our results for N — 1 and d ap- 
proaching 2 substantially deviate from the exact value 




FIG. 1: Illustration of the dichotomy procedure employed to 
identify the critical point S = S c . The present exemplary plot 
corresponds to v = 0.1, N = 1, and d = 2.6. For S > S c the 
plotted flowing coupling U/^'(p = 0) attains positive values 
at a finite scale, indicating flow to the symmetric phase. For 
S < S c the coupling u^ 1 {p — 0) decreases upon reducing the 
scale A as the system flows to the symmetry-broken phase. 
Fixed-point behavior represented by the plateau corresponds 
to the critical value <5 = <5 C a -0.01253855209. 
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FIG. 2: The correlation length exponent v(d, 1) plotted as 
function of system dimensionality d at iV = 1. The value 
u(3, 1) w 0.64 agrees well with the best available estimates. 
For d approaching 2 the obtained result substantially deviates 
from the exact value v exac t(2, 1) = 1. See the main text for a 
discussion of this discrepancy. 



X X 



X X 



x >< 



2/ e xact(2, 1) = 1. This discrepancy arises because the 
present exact RG truncation neglects renormalization of 
the momentum dependencies in the propagator. Indeed, 
the singularity of the propagator renormalization at crit- 
icality is quantified by the value of the anomalous dimen- 
sion r], which is relatively small in d = 3 (^3 « 0.036 for 
N = 1), but sizable (rj 2 = \ for N = 1) for d = 2. The 
present approximation, Eq. Q, obviously yields 77 = 0. 
This also perturbs the very structure of Eq. (fTU|) at d = 2, 
where the vanishing factor (2 — d) should be replaced by 
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FIG. 3: The correlation length exponent v(d, 2) plotted as 
function of system dimensionality d at N = 2. The value 
v(3, 2) « 0.68 compares well to the best available estimates. 
We observe a rapid growth of v(d, 2) upon reducing d toward 
2. This tendency is in line with the expectation that the XY 
model in d — 2 features an essential singularity of the correla- 
tion length when approaching criticality from the symmetric 
phase. 



(2 — d — rj) . For this reason we did not achieve true fixed 
point behavior at d = 2 and do not include the value 
of v precisely for d = 2 (see Fig. 2). This however does 
not create any obstacle in computing the Casimir forces 
in d = 2 and N = 1 close to the critical point in the 
symmetric phase (see the next section). 

For the same reason we are not able to resolve the 
critical behavior of the Kosterlitz-Thouless universality 
class (d = 2, N = 2), where nonzero rj is crucial to 
capture the physics even at a qualitative level. Inter- 
estingly, for N — 2, we observe a rapid increase of v(d, 2) 
upon reducing d towards 2. This may indicate the diver- 
gence of v(d, 2) and the onset of the Kosterlitz-Thouless 
physics, where the correlation length exhibits an essen- 
tial singularity in the symmetric phase. We also note, 
that the computed dependence v(d, 2) nicely fits a power 
law v(d,2) ~ (d - 2)~ A in the interval d e [2.15,2.4]. 
We obtain A « 1.31. We also note, that the values of 
v(d, N) in low d significantly improve upon including r/, 
i.e. going to higher order in derivative expansion. For a 
calculation resolving the (i-dependence for N — 1 at the 
next-to leading order in derivative expansion see Ref. HH. 



IV. CASIMIR FORCES 

Having resolved the critical behavior for the infinite 
system, we now consider finite L and proceed by com- 
puting the Casimir forces. The system state is chosen in 
the symmetric phase (S > S c ), but close enough to the 
critical point, so that the correlation length £ > 10 5 A ( ^ 1 . 
The L-dependent contribution to the free energy result- 
ing in the appearance of the Casimir force is identified 



<t(L) = lim 



F-fV 
D d-i 



where 



co y 

This leads to 

a(L) = k B TL [Etf_> (0) - U^ ( 



(13) 



(14) 



(15) 



The additional superscript in the effective potential 
Uj((p) indicates the finite size dependence that we set 
out to compute now. The potential U is evaluated at 
p = 0, because, as indicated above, we perform the com- 
putations in the symmetric phase. Utilizing Eq. ([5} with 
the cutoff given in Eq. ([7]) , we have 



dAU?( P ) = A d Z A d+1 
and 

d A u£( P ) = A d _ 1 L- 1 Z Q A d 



1 



N -1 



M-(p) M™(p) 



(16) 




(17) 



where n max — ^ and the quantities M ' (p, q 2 ) are given 
by Eq. © with the appropriate derivatives of the in- 
dependent potentials Ufc(p). All the dependencies on q 
dropped out after performing the momentum trace. We 
calculate the free energy a(L) given by Eq. (fTS"]) sup- 
plemented by the flow Eqs (fT6|) and (fT?]) , which are in- 
tegrated for a sequence of L-values. Eq. (fl5|) indicates 
that the thermal energy k B T sets the scale for a(L) 
and therefore also for the Casimir force. We plot the 
difference [J7^ (0) — U^ o (0)] versus L in the double- 
logarithmic scale, as exemplified in Fig. 4 for d = 2 and 
N = 1. The points tend to lie on a straight line for suf- 
ficiently large L, indicating scaling behavior. By fitting 
a straight line in the regime ^ > L > Aq" 1 , we find the 
slope -a = -2.0011 ± 0.0092 w -2 and the free coeffi- 
cient b = — 1.323±0.045. This translates to the computed 
Casimir force decaying as L~ a . The critical Casimir force 
amplitude Af is defined via 

ML) = ~ ~ -k B TA f L- a for large L , (18) 

and the general expectation is that a = d. The Casimir 
amplitude can be computed from the fit coefficient b. 
In the present case d — 2, N = 1 we find A/(2,l) = 
0.266 ± 0.012, which compares well to the exact result 

A exact = JL « 0.262. 

The above discussion exemplifies the procedure for the 
case d = 2, N — 1, where all results can be compared 
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FIG. 4: Logarithmic plot of the difference 

[l/iU Q (0) - tfJwo(O)] versus L for d = 2, JV = 1. 
Fitting a straight line to the scaling region yields the 
slope -o = -2.0011 ± 0.0092 and the free coefficient 
b = —1.323 ± 0.045. This translates to the computed 
Casimir force decaying as L~ a and the Casimir amplitude 
A f (2,1) = 0.266 ± 0.012, which compares well to the 
exact result A} xact = ^ ~ 0.262.^ Deviation from scaling 
behavior occurs for large and small L, where the condition 
^ > L > Aq" 1 becomes violated. 



to exact values. In principle, the present approximation 
is least reliable for d = 2, but nonetheless our result 
turns out to agree with the exact one. The method is 
obviously not restricted to this case and below we present 
our results for the Casimir force amplitudes A/(d, N) for 
N = 1 and N = 2 and varying dimensionality d between 
2 and 3. 



A. Ising universality class 

In this subsection we discuss the Ising universality class 
fixing N = 1, and perform the procedure outlined above 
varying d between 2 and 3. For these values of d, the ob- 
tained values of a which govern the decay of the Casimir 
force agree with the expected behavior a = d. The 
Casimir force amplitude as a function of d is plotted in 
Fig. 5. As already mentioned, the computed value of A/ 
in d = 2 agrees with the exact result A^ act = j^M- For 
d = 3 we obtain A f (3, 1) = 0.361 ± 0.026. This may be 
compared to numerical MC simulations of Ref. l35l . which 
give Af w 0.306. 

The variation of A.f(d, 1) is much more pronounced 
within the range d £ [2, 2.6] than for d > 2.6, where a 
constant value is compatible with our results within the 
errorbars. The relative error resulting from the fitting 
procedure outlined in the previous section is about 4% 
around d = 2 and increases to 7% around d = 3. The 
accuracy decreases for growing d because the quantity 
<t(L) computed via Eq. ([T5|) becomes smaller (at fixed 
L) and therefore the maximal value of L that can be 



FIG. 5: The Casimir force amplitude A/(d, 1) as function of 
dimensionality d in the Ising universality class N — 1. The 
extra point at d — 2 (square) corresponds to the exact result, 
and in d = 3 (triangle) to the prediction of MC simulations 
of Ref. |H. 

reliably handled numerically is lowered. 

Let us also compare the present result in d = 3 with 
previous RG studies relying on the e-expansion. The 
two-loop result of Ref. [25| retains terms to order e and 
yields Af (3, 1) « 0.22. As discussed in Ref. [26| (see also 
Ref. H3), the e-expansion is however problematic in the 
present case of periodic boundary conditions, yielding 
at three-loop level terms nonanalytic in e. The next- 
to-leading contribution is of the order e 3 / 2 and leads to 
Af (3, 1) corrected to w 0.39 in d = 3. 



B. XY universality class 

The computation proceeds along the same lines in the 
present case, where we put N = 2. This yields the extra 
contribution to the effective potential flow described by 
Eqs (HHJ) and (fTT)) resulting from the presence of the Gold- 
stone mode. As already discussed, the approximation 
level is too low to access the Kosterlitz-Thouless physics 
in d — 2 and we limit the range of d to d £ [2.2, 3]. The 
Casimir force amplitude Af(d, 2) is plotted as a function 
of d in Fig. 6. Similarly to the case N = 1, we observe a 
variation in the dependence A f (d, 2) for d < 2.6, while for 
d > 2.6 our results are consistent with A/(d, 2) = const. 
We find A/ (3, 2) = 0.707 ± 0.058, which may be com- 
pared to the MC valued Af c w 0.599. Similarly to the 
case N = 1, the MC result is lower and outside our er- 
rorbars. Our prediction is again closer to the MC results 
than the earlier perturbative RG studies, which yiel d 26 i 27 
Af (3, 2) f» 0.84. The dependencies depicted in Figs 5, 6 
are qualitatively similar, in particular we do not observe 
anything dramatic happening upon reducing d towards 
d = 2 in the 0(2) case (Fig. 6). 

It is also worthwhile noting the relatively large error- 
bars resulting from the fitting procedure illustrated in 
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FIG. 6: The Casimir force amplitude A/(d, 2) as function of 
dimensionality d in the XY universality class N — 2. The 
extra point in d — 3 (triangle) corresponds to the prediction 
of MC simulations. 



detail at the beginning of the section. The errors account 
only for the numerical inaccuracies. We have nothing to 
say here about the errors caused by the approximation 
made while deriving the flow equations. As already re- 
marked, the errorbars are wider for d &s 3 than in d ~ 2. 
The inaccuracies are larger than one might perhaps ex- 
pect, because calculation of the Casimir amplitude re- 
quires finding the free coefficient in the fit straight line 
given only points far from log(LAo) = 0. For this reason 
a small tilt variation causes sizable changes in the point 
of interest, where the fit line intersects the y-axis. This 
results in significantly worse accuracy in determining the 
free coefficient than the slope. 



V. SUMMARY AND DISCUSSION 

In this work we have performed a renormalization- 
group study of Casimir forces occurring in the classical 
rf-dimensional O(N) models confined to a slab- like geom- 
etry and subject to periodic boundary conditions. In this 
case the system preserves translational symmetry and the 
Casimir forces arise solely due to fluctuations. These are 
long-ranged and universal at criticality. Our approach re- 
lies on a truncation of the one-particle- irreducible variant 
of functional RG at leading order in the so-called deriva- 
tive expansion. We believe this approach complements 
the MC calculations in that the accessible system sizes 
are by far larger and no reference to finite-size scaling 
theory is required. On the other hand, it does not suffer 
from complications occurring in perturbative RG treat- 
ments and related to nonanalyticities in the e-expansion, 
revealed beyond two loops. It is also applicable both in 



d close to 2 and 3, and we scanned a range of d- values in 
our calculations. 

Our results agree with the expectation that the criti- 
cal Casimir force decays with the separation L between 
the walls as L~ d and the thermal energy fcgT sets the 
relevant scale. We have computed the critical Casimir 
force amplitudes A/(d, N) for the Ising (N = 1) and XY 
(iV = 2) universality classes varying dimensionality be- 
tween 2 and 3 and remaining in the symmetric phase, 
but very close to the critical point. Our findings in d = 2 
compare well to the exact results. In d = 3 the com- 
puted Casimir amplitudes reasonably agree with numer- 
ical simulations, they are also closer to MC than the ear- 
lier field-theoretic RG calculations. The level of agree- 
ment (in particular in d = 2 and N — 1) suggests that 
the details of the momentum dependencies of the corre- 
lation function are not relevant for the magnitude of the 
Casimir forces at criticality. This is somewhat surprising, 
considering that the presently studied effect arises exclu- 
sively due to critical fluctuations. Indeed, our calculation 
neglects the nonzero anomalous dimension rj whatsoever 
and also gives incorrect values of the correlation length 
exponent v in d = 2. Still, the results agree well with the 
exact findings in d = 2. 

From a principal point of view, the present approach 
is applicable for all the O(N) models at arbitrary di- 
mensionality d with periodic boundary conditions. One 
exception is the Kosterlitz-Thouless case, which requires 
going to higher order in derivative expansion. In practice 
however, the computation of the Casimir forces requires 
increasing numerical accuracy with growing d, and may 
become problematic already at d = 4. Another interest- 
ing regime is N — > 00, where one may again compare to 
exact results. It is however expected, that the Casimir 
amplitudes become small in this limit and computations 
proceeding along the present line might yield unreason- 
ably large relative errors. 

It might be interesting to extend the present work 
to account for boundary conditions explicitly breaking 
translational symmetry. This problem is of interest also 
from the point of view of comparison to experiments per- 
formed on helium and fluid mixtures. The other inter- 
esting extensions which should be within the range of 
the present approach include computations within the 
ordered phase, where the soft fluctuations arise due to 
continuous symmetry breaking, and extracting the scal- 
ing functions for the Casimir force. 
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